Retarded radiation from colliding black holes in the close limit 
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We use null hypersurface techniques in a new approach to calculate the retarded waveform hom 
a binary black hole merger in the close approximation. The process of removing ingoing radiation 
from the system leads to two notable features in the shape of the close approximation waveform 
for a head-on collision of black holes: (i) an initial quasinormal ringup and (ii) weak sensitivity to 
the parameter controlling the collision velocity. Feature (ii) is unexpected and has the potential 
importance of enabling the design of an efficient template for extracting the gravitational wave 
signal from the noise. 
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I. INTRODUCTION 

We compute the retarded radiation produced by the 
collision of two black holes in the close approximation 
where the collision occurs in the distant past. The 
retarded wave corresponds to the perturbation of a 
Schwarzschild spacetime in which there is no ingoing ra- 
diation. The calculation is carried out by means of a 
characteristic evolution algorithm. In previous work § 
we applied characteristic evolution to determine the close 
approximation retarded waveform from a white hole fis- 
sion. Here we find that the retarded waveform from a 
black hole merger has distinct features which do not ap- 
pear in the fission waveform: (i) an initial quasinormal 
ringup and (ii) weak sensitivity to the parameter con- 
trolling the collision timescale. Feature (ii) has potential 
importance for the design of an efficient template for ex- 
tracting the gravitational wave signal from noise. 

The close approximation |^|, based on the perturba- 
tion of a stationary single black hole spacetime, has been 
extremely useful for understanding the binary black hole 
problem by means of Cauchy evolution |, |, |, |, |, 
y, |l^, . Our previous application of characteristic 
evolution to compute the retarded solution for a white 
hole fission [Q is the starting point for the present work. 
Under time reversal, the retarded solution for a white 
hole fission is equivalent to an advanced solution for a 
black hole merger. In the advanced solution, ingoing ra- 
diation is absorbed by the black holes but no outgoing 
radiation is emitted. The advanced solution Q provides 
stage I of a two stage approach to the binary black hole 
problem |jl^, ^ |T^. In stage II, treated here, we present 
the physically appropriate retarded waveform for a black 
hole collision corresponding to the retarded solution in 
which there is no ingoing radiation from past null infin- 
ity I~ in the background Schwarzschild spacetime. 

These perturbative results furnish a reference point for 
the physical understanding of a fully nonlinear version of 
the same approach [|l3[ ^ which is being pursued us- 
ing existing characteristic evolution codes |l^ . These 
results also complement the physical picture obtained in 



this approximation by Cauchy evolution. 

The advanced solution for a black hole merger is sim- 
pler to compute in a characteristic approach because of 
the relative location of the two null hypersurfaces on 
which boundary information is known. One of the null 
hypersurfaces is the black hole event horizon 7^+, whose 
perturbation corresponds to the close approximation of 
a binary black hole. In the advanced problem, the other 
null hypersurface (in a conformally compactified descrip- 
tion) is future null infinity I"*" where the outgoing ra- 
diation vanishes. Since TL^ and I"*" formally intersect 
at future time infinity they can be used to pose 
a double-null initial value problem jl^, |20[ ^ to 
evolve backward in time and compute the exterior region 
of spacetime. (Potential difficulties in dealing with /+ 
are avoided by posing the problem on an ingoing null hy- 
persurface i/"*" which approximates X"*" and intersects Ti^ 
at a late time when the perturbation has effectively died 
out.) On the other hand, in the retarded problem, the 
other null hypersurface is past null infinity I~ where the 
ingoing radiation is required to vanish. Because Ti/^ and 
T~ are disjoint, there is no direct way to base a charac- 
teristic initial value problem on these two hypersurfaces. 

The advanced solution provides the ingoing radiation 
incident from X" . In the context of linear perturbation 
theory, the time translation and refiection symmetries of 
the background provide a mathematically well defined 
way to carry out stage II by using this ingoing radia- 
tion to generate a "source free" advanced minus retarded 
solution. A retarded solution is then produced by su- 
perposition with the stage I solution. In this paper, we 
deal with the subtleties produced in implementing this 
procedure computationally. 

Time reflection symmetry allows the merger of two 
black holes to be described as the fission of a white hole. 
It is convenient here for computational purposes to for- 
mulate the binary black hole problem in this time re- 
versed sense as a white hole fission since the characteristic 
evolution then takes the standard form of being carried 
out forward in retarded time. From a time reversed point 
of view, the stage I solution is equivalent to the retarded 
solution for a white hole fission, with the physically ap- 
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propriate boundary conditions that radiation is emitted 
but that there is no ingoing radiation from I~; and, un- 
der time reversal, the stage II solution is equivalent to 
the advanced solution for a white hole fission, with the 
boundary condition that radiation from X~ is absorbed 
but no outgoing radiation is emitted. 

In Sec. II, we review the background material upon 
which this work is based 0|. We use a characteristic 
reformulation of the Teukolsky equations |^ governing 
the perturbative Weyl tensor. Previous work dealt with a 
characteristic evolution algorithm for solving the Teukol- 
sky equation for the ingoing null component of the Weyl 
tensor |Q| . In order to carry out Stage II, we extend this 
work to an evolution algorithm for evolving the outgo- 
ing null component. Characteristic data obtained for the 
nonlinear description of a binary black hole spacetime is 
used to induce close approximation data for the pertur- 
bative solution considered here. 
In Sec. 
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we provide the technical details of our 
global strategy for producing the retarded waveform from 
a black hole process. In Sec. IV, we describe the numer- 



ical implementation and testing of the evolution algo- 
rithm. Because of difficult computational problems of 
a global nature, we develop two independent evolution 
algorithms and carry out extensive testing in order to 
demonstrate that the final waveform is free of numerical 
artifacts. 

In Sec. we present the properties of the retarded 
waveform from a head-on black hole collision in the close 
approximation. It has been suggested that the black hole 
merger waveform might have a distinct signature that is 
roughly independent of the details of the initial condi- 
tions. Our results do confirm the generally accepted view 
that the waveform has a quasinormal ringdown stage 
characteristic of the final black hole mass (which is here 
the background mass). However, we find that the fea- 
tures of the main waveform are sensitive to the pres- 
ence of ingoing radiation fields of the type that would 
be present, for instance, in the evolution of time refiec- 
tion symmetric Cauchy data. 

Notation and Conventions: We use a metric of sig- 
nature ( V H — h) and a null tetrad with normalization 

l°-na = -m'^fha = -1, so that gab = 2{m(^amb) - l(anb))- 
We introduce standard angular coordinates = (6*, 0) 
to represent the unit sphere metric qAB and set q^^ = 
q(^q^\ where q^^qsc = "^i^h q^ = (l,z/sin6'). We 
use q"^ to define the 3 operator with the convention 9/ = 
q^dAf, for a spin-weight function /. Complex conju- 
gation is denoted with a "bar", e.g. 3?(/) ^ {f + f)/2. 
These conventions lead to a different form of the pertur- 
bation equations from that originally given by Teukolsky. 



II. BACKGROUND FORMALISM 

Various coordinates were used in Ref. to describe 
curvature perturbations of a Schwarzschild background. 
The simplest global choice are Israel coordinates ||2^ in 



which the Schwarzschild metric takes the form 

ds2 ^ ^Wdu'^ - 2dudX + r'^qABdx^dx^, (2.1) 

where u is an affine parameter along the ingoing null 
hypersurface r = 2M that forms the white hole horizon, 
A is an affine parameter along the outgoing null geodesies 
in the direction of increasing r — 2M — \u/{AM) and 
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2A2 



(2.2) 



(Affine freedom has been used to set A = on the white 
hole horizon, to set u = at the r = 2M bifurcation 
sphere between the black hole and white hole horizons 
and to set g"'^ {dau)di,\ = —1). Israel coordinates cover 
the entire Kruskal manifold with nonsingular and simple 
analytic behavior. The exterior quadrant is given by u < 
0, A > 0. 

Other coordinates more suitable for certain numerical 
purposes are outgoing Eddington-Finkelstein coordinates 
(m, r), in which the Schwarzschild metric is 



ds^ = - 1 )du 

r 



,~,2 



2dudr + r^dn'^. 



(2.3) 



Here u = —Me and u is a Bondi retarded time co- 
ordinate related to the standard Schwarzschild time co- 



ordinate t hy u = t 



with r* 



2Mlog( 



2M 
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These outgoing Eddington-Finkelstein coordinates patch 
the lower half of the Kruskal manifold composed of the 
exterior quadrant and the quadrant above the initial sin- 
gularity, with u oo at the black hole horizon. 

Israel coordinates induce the nonsingular complex null 
tetrad 



r 



'V^'u^i — r^ [0,1,0,0], 
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The Teukolsky equations |22[| for the vacuum perturba- 
tion of the outgoing null component of Weyl curvature, 
ipo — Cabcdl'^'m''l'^m'^ (in Newman-Penrose notation [p4|) 
and the ingoing null component 7/^4 = Cabcdn°"fn^n'^'iTi'^ 
are 

1^2 j^2 

(^0 + ^7^)^0 = 0, and + = 0, (2.6) 
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where 
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is the Laplacian defined by the metric ds^ = — (1 — 
2M/r)du^ — 2dudr induced by the backgrou nd o n the 
2-dimensional (w, r) subspace. Equations ( 2.10| ) and 
( ^.ll|) asymptote to one-dimensional wave equations for 
solutions whose radial derivative falls off uniformly as 
0{l/r'^). The limit of F4 determines the outgoing grav- 
itational radiation waveform while the limit of _fo is re- 
lated to the retarded quadrupole moment of the system. 
More precisely, limr'04 = duN, where N is the standard 
definition of the Bondi news function. 



J r2 ;^?s • ii 1 i J In these variables, the deviation of Eq's. (2.8) and 

and L = —00 IS the angular momentum squared opera- ij-™„-Ji „ i-„ - •„ jT, „j JT 

tor. 

The Weyl components ipo and ip4 have spin-weight 
s — 2 and s — —2, respectively. It is useful to convert 



Eq's. (2.6) - (2.7) into spin-weight-zero equations by set- 
ting tpo — d^^Q and tp4 = §^$4. Furthermore, the com- 
putation of asymptotic behavior is enhanced by working 
with the fields Fq — r^$o and F4 = r<l>4, which gener- 
ically have finite non-vanishing limits at X+ for asymp- 
totically flat perturbations. 

Although the Weyl scalars formed from the (/",«") 
null vectors are non-singular in the full Kruskal mani- 
fold, they are not well adapted for numerical accuracy. In 
order to compute radiation at late times near 2""*" it is ad- 
vantageous to consider a boosted tetrad (l°-,h°-, m°, to°), 
with 1°- ~ — and satisfying I'^ha — —1, which is 
adapted to Bondi time u at X"*" rather than the affine hori- 
zon time u. We denote the corresponding Weyl compo- 
nents by ■00 = Cabcdl°"rn''l''m'^ and Tp4 = Cabcdn'^'fnf'n'^'fn'^ ■ 
In particular, the late time behavior of the radiation 
waveform can be more accurately computed by an evo- 
lution algorithm for the Weyl component ip^. 

Setting ■ipQ = S2|>o and xpi = the Teukolsky 

equations for the fields Fq = r^$o and F4 — r$4 reduce 
to the spin-weight zero forms 

{D^ + So)Fo = 0, (2.8) 
{D^ + S4)F4 = 0, (2.9) 



where 



and 



16M(r-3Af) 30Af (6-1-99) 

JO — n dx n \ 5 

= _4(. 3A/)^^_30M^(6 + 99) ^^^^^^ 



16Af(r - 3Af) 6Af (2 + 99) 
^ "A ^ H ^ — 
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from a 1-dimensional wave equation is independent of 
time at a fixed r. Because F4 = (u/4Af)2F4, where F4 
is regular throughout the Kruskal manifold (since it is 
constructed with a regular basis), it follows that F4 ^ 
as the black hole horizon is approached. This facilitates 
an accurate long term evolution of the waveform using 
Eq. ( pj| ). 

In contrast, Fq = {4:M/u)'^Fo, so that Fq is singular 
on the black hole horizon, and thus a poor choice of vari- 
able for long term evolution. The opposite signs of the 
coefficients of dr in Sq and S4 are responsible for this 
behavior, as can be seen by ignoring the remaining po- 
tential terms and freezi ng t he coefficient of dr a,t r — 2M, 
so that Eqs. (2.8) and (2.9) reduce to 



{2du - dr - ^)drFo = 0, 

{2du-dr + ^)drF4 = 0, 



(2.13) 
(2.14) 



in terms of retarded Bondi coordinates. In this approx- 
imation, both of these equations admit purely outgoing 
waves F(u). However, an ingoing _Fo wave has the expo- 
nentially singular behavior 



#0 = f{u + 2r)e- 



(2.15) 



as an initial pulse /(mq + 2r) approaches the black hole 
horizon as u — > 00. An ingoing F4 wave behaves in op- 
posite fashion and decays exponentially on approach to 
the black hole. 

and (12. 9|) make it clear 



The different forms of Eqs. (2 



that the Weyl component V'o in the outgoing null direc- 
tion ?° and the Weyl component "04 in the ingoing null di- 
rection n°' are not related in a way which makes manifest 
the Schwarzschild time reflection symmetry T, defined by 
T{t, r, 9, (p) = {—t, r, 6, (/>) in Schwarzschild coordinates or 
by T(u, r, 9, (f>) — (— w, r, 9, (f>) in terms of Bondi retarded 
and advanced times u = t ~ r* and v = t + r* . The 
time reficction symmetry could be incorporated explic- 
itly into the pair of Teukolsky equations by introducing 
null tetrad vectors = and N'^ — {l/a)n°' satis- 
fying TL° = —N'^. However, the explicit form of the 
required boost, 
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8M^ - Xu 
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makes it clear that such a time symmetric formulation 
would introduce singular behavior at both the black and 
white hole horizons. 

This time symmetric tetrad is useful for formulat- 
ing the time reflection properties of solutions of the 
Teukolsky equations using other tetrads. Let ^'4 = 
CabcdN°'mJ'N'^rh'^ = '^{u,r,6,(j)) be a perturbative so- 
lution for ^'4. Then the time reflection symmetry im- 
plies that *o = CabcdL^m^L'^m'^ = ^{-v,r,9,4>) is a 
perturbative solution for ^q. This correspondence maps 
a retarded solution (no incoming radiation) for \E'4 into 
an advanced solution (no outgoing radiation) for vPq. In 
terms of the 1°" and n° Weyl components, the solution 
^^4 = tlj{u,r,9,(j)) corresponds under time reflection to 
the solution 



■00 



4r^ 



(r - 2M)' 



(2.17) 



A. Close approximation data 

Data for the double-null formulation of the characteris- 
tic initial value problem are given on a pair of intersecting 
null hypersurfaces. The data for a fissioning white hole is 
posed on the ingoing null hypersurface 7i~ and an early 
outgoing null hypersurface which extends from 7i~ 
to X+. In Ref. iQ, consistent double- null data for ip^ 
for a fissioning white hole was obtained from the metric 
data for the nonlinear version of the problem in terms 
of a spin-weight-2 field J describing the inner conformal 
geometries of Ti~ and Q. 

The underlying field J for a white hole fission is pro- 
vided by a conformal model based upon an ingoing null 
hypersurface Af emanating from a prolate spheroid em- 
bedded in a flat space j2|, |2^. Let (f,6,(f>) be stan- 
dard spherical coordinates for the inertial time slices 
t = constant of Minkowski space. When the eccentricity 
of the spheroid vanishes, the Minkowski null hypersur- 
face reduces to the light cone from a sphere t — f — a. 
The intrinsic geometry of M is described by a degenerate 
rank-2 metric r^hAB where J = q^q^hAB ■ To linear 
order in the eccentricity, the conformal model gives 



jii,d) 



t — a 



(2.18) 



This conformal geometry of JV models the conformal 
geometry of the white hole horizon in the close limit 
where the fission occurs in the infinite future. The area 
coordinate r of the geometry oiTi^ differs from the cor- 
responding f of the geometry of J\f in order to yield an 
asymptotically constant surface area for the initial white 
hole in the infinite past. For vanishing eccentricity, this 
produces a Schwarzschild horizon and a small eccentricity 
generates the close approximation. 

The Raychaudhuri equation, which governs the expan- 
sion of a null hypersurface, forces a difference in the afhne 
parameters on Af and 'H~ . With a suitable choice of 



affine scale, the relation between i and the afhne param- 
eter u on the white hole horizon is 



d^I - = (3 



f2(f 



5f + f2)2 



(5 - V13) - 2f 
(5 + yiS) - 2t 



(2.19) 



where t ~ i — a. At early times Eq. ( p. 19 ) implies u ^ i 
but as the Minkowski null cone pinches off at t = a 
the corresponding afhne time on the white hole horizon 
asymptotes to w ^ cx). In terms of the inverted pair-of- 
pants picture for a white hole fission, the pants legs are 
mapped to m = 00, so that in the close limit the individ- 
ual white holes are mapped to future infinity along the 
white hole horizon in the Kruskal manifold. The details 
are discussed elsewhere in a treatment of fully nonlinear 
null data for the general two black hole problem ^ . 

Close limit data for J(u, 6) on the w hite hole horizon 
is determined by integrating Eq. ( 2.1S| ) and substituting 
into Eq. ( 2.18 ). In order to eliminate nonessential inte- 
gration constants, we set u = at the r = 2M bifurcation 
sphere. Then, up to scale, the close data are determined 
by the single parameter 



V = 



J\u=0 
'^J\u—oc 



T\u=0 



(2.20) 



where 77 > is a scale invariant parameter determining 
the yield of the white hole fission. (In the time reversed 
scenario of a black hole collision, r/ is related to the in- 
elasticity of the collision.) In the linear approximation. 



V'4 



2 J .uu • 



(2.21) 



In this way, the conformal horizon model determines 
the close limit data for ip4 on the horizon 7i~. By choos- 
ing J'~ to be an early outgoing null hypersurface approx- 
imating X~ , we prescribe data ■04 = on JT" , which is 
consistent with the asymptotic falloff -04 = 0{u~^) along 
Ti^ implied by Eq. (2.18). This approximates the condi- 
tion that there be no ingoing radiation from X" . 

In terms of the spin- weight-zero potential -^4, the close 
data on Ti.^ have the normalized form 



-iAdf){Adf)- 



(2.22) 



after factoring out the £ — 2 angular dependence. 

The time dependence of the close approximation data 
is quite mild when expressed as a functio n of f, as 
in Eq. ( 2.22 ). However, the relationship ( 2.19 ) pro- 
duces a sharp transition region where the behavior of 
f{u) changes from the asymptotic form dr/du — > 1 as 
f —00 to dr/du ^ as f ^ Q. For large values 
of the parameter i], this produces sharply pulse shaped 
data |l). 
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III. COMPUTING RETARDED WAVEFORMS 
FROM BLACK HOLES 



A. Removing the ingoing radiation 



data would be artificial. In stage II, the conversion from 
an advanced to retarded binary black hole solution au- 
tomatically introduces this ringdown in a natural way in 
addition to removing the ingoing radiation. 



In the time reversed scenario being pursued here, our 
strategy for computing the retarded solution for a black 
hole collision is based upon the computation of the ad- 
vanced solution from a white hole fission. Stage I of the 
approach, i.e. computation of the retarded solution for a 
white hole fission, has already been carried out . Close 
approximation data for ■04 provided by the conformal 
horizon model was evolved to determine the stage I so- 
lution ■04^ throughout the exterior region of spacetime. 
In practice, the evolution computes the corresponding 
spin-weight-0 potential F^j, whose value on X"*" deter- 
mines the outgoing radiation field F^j^out- In stage I 
of the white hole problem, there is no ingoing radiation 
from . In stage II, the strategy for converting the re- 
tarded solution to an advanced solution with no outgoing 
radiation is to superpose a perturbative solution whose 
outgoing field equals —Fij^ouT- 

This procedure for converting from an advanced to a 
retarded field in a white hole spacetime differs from the 
standard procedure in a flat spacetime. The retarded 
solution ^RET of the flat space wave equation = 
with source S may be converted into the advanced solu- 
tion ^ADV by superposing the source free (homogeneous) 
solution H ~ ^ADV — ^ret- In the vacuum white hole 
(or black hole) problem, there is no source in the exte- 
rior spacetime. The white hole horizon itself acts as the 
source of the exterior solution. In this vein, specification 
of the "source-free" solution which transforms a retarded 
white hole perturbation to an advanced white hole per- 
turbation requires specification of horizon data on Ti"*" as 
well as the outgoing data —F^^j^qut at X"*". This hori- 
zon data must not introduce any features which change 
the physical interpretation of the system as a white hole 
fission. For example. It is essential not to superpose a 
black hole merger on the white hole fission, as would re- 
sult in the specification of a Cauchy problem with time 
symmetric data. 

We carry out the conversion to an advanced solution 
by evolving backward in time the double null problem 
with data i^4 = on the black hole horizon 7^+ and data 
F4 = —Fij^ouT at X"*". This data leaves the black hole 
horizon unchanged and its evolution determines the ana- 
logue of the "homogeneous solution" H. The superposi- 
tion Fi ji — F41 + H then produces the desired advanced 
solution for a white hole fission with no outgoing radia- 
tion. The stage II solution leaves the black hole horizon 
unchanged but perturbs the white hole, essentially mak- 
ing it ring. In the time reversed black hole scenario, this 
is a physically desirable effect because the binary black 
hole data provided by the conformal horizon model has 
a power law decay in u, which corresponds to a purely 
exponential decay in u without quasinormal oscillations. 
A priori introduction of ringdown in the final black hole 



B. Producing the outgoing waveform 

The Weyl components ipo and 7/14 form the major part 
of the geometrical description of a Schwarzschild pertur- 
bation. The waveform at X"*" is described directly by "04 • 
Stage II provides indirect waveform information in terms 
of -00 J which can then be used to obtain the waveform 
by means of differential relations between tjjQ and "04 ■ 
Asymptotically, at X+ , these differential relations can be 
obtained either from the vacuum Bianchi identities 
or from the asymptotic behavior of the metric |2^, |2^, pO|| . 
Here we present the latter metric approach, since it pro- 
vides a complete asymptotic description of Schwarzschild 
perturbations. 

In the {u,r,x^) Bo ndi coordinates with background 
Schwarzschild metric (2.3), a vacuum perturbation near 
X+ perturbs the metric in the asymptotically flat 
form iotl 



Sgar- 


= SgABQ^q^ = 0, 


(3.1) 


Sguu 


= — + 0(^3)' 


(3.2) 


SguAq^ 


Be 2C ^, 1 , 
4 3r 


(3.3) 


SgAsq^q^ 


= Cr+^ + 0(l), 


(3.4) 



where M and L are spin- weighted functions of 

[ii^x^). The spin-weight-2 functions C and Q comprise 
the asymptotic part of the initial null data, in metric 
form, on a constant u hypersurface. In Bondi's termi- 
nology, the spin-weight-0 function M. is the (perturbed) 
mass aspect of the system and the spin- weight 1 function 
C is the angular momentum aspect. The time deriva- 
tive of SgABQ^q^ can be determined from the evolution 
equations. The time derivative duC is the Bondi news 
function and requires a radial integration to determine 
its value in terms of a boundary condition on an inner 
worldtube. The time derivative of Q can be determined 
locally. Similarly, the time derivatives of A4 and £ can 
be determined locally from the components of Einstein's 
equations representing conservation laws. In an asymp- 
totic Bondi frame, the relevant evolution-conservation 
equations are Eq's (7.19) - (7.22) of Ref. They imply 



duQ 

duM 
duC 



ac Mc 

3^2' 

g^c -1-82(7 
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9(3^(7 -g^c) 

16 ■ 



(3.5) 
(3.6) 
(3.7) 
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The asymptotic values of the Weyl components are ob- 
tained from the metric according to 



IV. NUMERICAL ALGORITHMS AND 
ACCURACY 



3 

1 



(3.8) 
(3.9) 



This allows determination of ■04 in terms of ipQ by sol ving 
for C in terms o f Q. By ta king 2 -i-derivatives of Eq. (3.5) 
and using Eqs. (3.6) and ( |3.7D to ehminate M. and C, we 
obtain 



dlQ = ^^B'd,C + ^diC. 



(3.10) 



Replacing C and Q by spin- w eight potentials defined 
by C = 9^c and Q — Eq. ( 3.10 ) reduces to 



(3.11) 



When c is real, the solutions of Eq. (3.11) for which 
q = are solutions to the linearized Robinson- Trautman 
equation. Since the Robinson- Trautman solutions are 
singular as u —00, the requirement of nonsingular 
initial conditions leads to a one-to-one correspondence 
between c and q. The axisymmetric head-on collision 
of two black holes corresponds to the case where c is a 
real {£ = 2, to = 0) spherical harmonic. In that case, 



Eq. ( p.llD gives 



duC 



M 



(3.12) 



When c is complex, the solutions for q = correspond 
to a twisting version of a Robinson- Trautman perturba- 
tion. In this case, the imaginary part of c goes expo- 
nentially to cxD as u — > -l-oo and the relation between c 
and q is more complicated. The complex case will not be 
considered further here. 

For the quadrupole case, Eq. (3.11) determines the 
asymptotic relation at X"*" 



1 



(3.13) 



which determines the waveform tj^l = "d^F^ at X+ in 
terms of the asymptotic field -00 — S^Fq. This relation 
can also be expressed in the integral form 



^5 = duN 



(3.14) 



1 1 - 

= 5^^-'"/'' /" die^'/'"dfm, (3.15) 



where N is the Bondi news function. 



The stability of the background Schwarzschild geom- 
etry and the insensitivity of the Weyl tensor to gauge 
problems imply the absence of any theoretical limita- 
tions to the accurate numerical evolution of the Teukol- 
sky equation in the exterior Kruskal quadrant, where the 
background spacetime is nonsingular. However, there is 
a practical limit to how close a numerical evolution can 
produce meaningful results near future time infinity 
or past time infinity /~ . There is no nonsingular way to 
treat time infinity as a boundary in a conformally com- 
pactified manifold, as is possible for null infinity. Thus 
time infinity can only be approached asymptotically and 
the asymptotic values of fields can depend on the direc- 
tion of approach. This direction dependent limit at 
makes an evolution of ■00 difficult. 

Numerical difficulties have led us to implement the 
Stage II evolution of Fq using two distinct approaches, 
one based upon a retarded time foliation and the other 
on an advanced time foliation. Both approaches require 
finding variables which allow for stable evolution of waves 
propagating inside the peak of the Schwarzschild poten- 
tial, which turns out to be difficult. The retarded time 
algorithm uses an evolution variable which is adapted to 
the falloff at X"*", and the advanced time algorithm uses 
an evolution variable adapted to the falloff at X~ . Both 
approaches have similar accuracy but the advanced time 
algorithm is more efficient at treating behavior near Ti"*" 
and the retarded time algorithm is more efficient near 
X+. 



A. Retarded time algorithm 

We base the retarded time algorithm on the evolution 
variable 



Fq = {l-2M/rfFo, 



(4.1) 



which is finite on 7i+ , as well as on 2'^ where the wave- 
form is read off. However, the factor (1 — 2M/r)^, which 
makes Fq vanish identically on 7i~, introduces a singu- 
larity at Ti.^ in the associated Teukolsky equation. In ad- 
dition, _Fo is infinite on X~. We deal with these problems 
by initializing the double-null problem on an ingoing null 
hypersurface IC~ approximating Ti.^ and on an outgoing 
null hypersurface approximating X^ . 

In (u, v) double null coordinates, the Teukolsky equa- 
tion for Fq {£^2) is 



3A'f (1 - 2M/r) 
2^ 



^^n = 0. 



(4.2) 



The code is implemented in coordinates {u, p), where p = 
arctan(^i;/a) and a is an adjustable parameter (typically 
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a = 40M). These coordinates compactify X+ at p = 7r/2. 
The Teukolsky equation takes the form 



dudpFo = - 



2 (r - 2M) 



3 a (r - 2M) (l + tan^ p) 



Fo, (4.3) 



which is nonsingular in the region — Tr/2 < p < tt/2. 
However, compactification of X"*" using double null coor- 
dinates necessarily leads to a loss of d ifferentiability at 
X+ The coefficients in Eq. (4.3) have finite lim- 

its at X+, but the OuFq coefficient is not differentiable 
at 2"+. In {u,l = l/r) coordinates Fq can have analytic 
dependence on / at but the transformation to (w, p) 
coordinates introduces logarithmic behavior so that Fq 
is not twice differentiable with respect to p at X+. Our 
tests show this still allows the computation of Fq on X+ 
to 2nd order accuracy in grid size. Similar logarithmic 
behavior exists at I"*" in our original Stage I algorithm 



We implement Eq. (4.3) as the first differential order 
system 



, , 2(r- 2M) - 
9fiFo = g{u,p)-^ '-Fa 



(4.4) 



dp9 



2a (l -I- tan^ p) 



a (l + tan^ p) 



2r4 

X (r-2Af)(3M-h2r)Fo, 



(4.5) 



where g is introduced as an auxiliary variable. 

The numerical [u, p) grid consists of points on the u = 
const foliation. For reasons of economy and accuracy, 
the grid points are uniformly spaced in v in the region of 
compact support where the initial data is non-zero and 
are uniformly spaced in p in the out er r egion where the 
initial data vanishes. We solve Eqs. (4.4) and (4.5) using 
a 2nd order accurate Runge-Kutta scheme. 



B. Advanced time algorithm 

We base the advanced time algorithm on the evolution 
variable 

Fo = (1 - 2M/r)2r$o = ^i^o, (4.6) 

which is finite on and on ?i+ . Since Fq automatically 
vanishes on 7i~, we must again initialize on an ingoing 
null hypersurface K,~ approximating 7i~. In addition, Fq 
vanishes on I"*" . Consequently, the signal on X+ must be 
obtained by transforming Fq to Fq on successively far- 
ther out V = const slices and extrapolating to v = oo. 



Although the advanced time algorithm can be initialized 
on X~ this gives no distinct advantage in computing the 
waveform at X+ and, in order to facilitate comparison 
with the retarded time algorithm, we initialize on an out- 
going null hypersurface J^^ which approximates . 
The advanced time evolution variable Fq satisfies 



1 n _L 2(r-3M) 

OuOy + du 



(1 - 2M/r)(2r + 3Af) 
2^ 



Fq^O. (4.7) 



This is precisely the equation for F^ under the transfor- 
mation V ^ —u. (Under time reflection, Fq — ATF4). 
Thus the advanced time algorithm is equivalent to an 
algorithm for evolving F^ backward in time. 

The evolution proceeds along v = const null hypersur- 
faces, which are compactified by the coordinate transfor- 
mation u ~ atan/i, where — 7r/2 < /i < 7r/2 in the range 
from to 7i+, and a is again an adjustable parame- 
ter (typically a = 200 A/ in the advanced time case). In 
these coordinates, Eq. (4.7) becomes 



a(l-2M/r)(2r + 3Af) 
2r^ cos^ p. 



Fq^O. (4.8) 



This equation is non-singular in the region —Tr/2 < p < 
tt/2, but the coefficient of Fq is not differentiable on 
[p = — 7r/2). On an ingoing null hypersurface, Fq can 
be expressed as a power series in I = l/r. This implies 
that the associated series in p contains logarithmic terms 
of the form /i(l + plogp), and is not twice differentiable 
(with respect to p) on X~ . Our tests indicate that the 
code is 2nd order convergent despite this mild lack of 
differentiability. 

We implement Eq. (4. 
form 



in the first differential order 



d^,FQ 
dv9 



2{r-iM)g 

a{r - 2M){2r + 3M)Fq 
2r^ cos^ p ' 



(4.9) 



(4.10) 



where we have introduced the auxiliary variable g. The 
numerical {p, v) grid consists of points on the v — const 
folia tion equidistant in p. The hypersurface equation 
(Oh is solved by 2nd order accurate integration and Eq. 
( f.lOj) is solved by a 2nd order accurate Runge-Kutta 



scheme. 



C. Extracting the Retarded Waveform 

The retarded waveform F4 on I"*" is obtained from Fq 
by means of Eq. ( 3.131 ). Given the values of Fq, we 
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integrate Eq. (3.12) by the midpoint rule to obtain F4, 
using the initial condition that -F4 = 0. 

The retarded algorithm provides Fq on X+ directly. In 
the advanced algorithm, where X"*" is not in the evolu- 
tion domain, we extrapolate the waveform at X"*" from 
a quadratic polynomial based up on ap proximate values 
of F4 obtained by applying Eq. (3.13) on three ingoing 
null hypersurfaces approximating X+ (where we substi- 
tute r'^Fo for Fq). The three null hypersurfaces consist 



of the advanced time slices Vf, 
(typically) if = 4000M. 



Vf/2, and v-f/3, where 



D. Initializing the Stage II evolution 

Initial data for the stage II evolution, i.e. the value of 
ipo on , is obtained by a time reflection of the retarded 
Stage I solution for ip4 on Technical problems arise 
here because the stage I evolution does not provide ac- 
curate data at late times on J'~^ Q. We deal with this 
by smoothing the Stage II initial data to zero outside a 
region of compact support where the data would other- 
wise be small. This modification of the initial data also 
guarantees consistency with the boundary data ipo — 
on /C~, which approximates the boundary condition that 
the perturbation vanishes on 7i~. 

duN onX+ requires that the integral 



The relation V'? 



(4.11) 



vanish since the news function N must vanish at either 
limit in order that the radiated mass loss be finite. By 
time reflection, it follows that the integral 



dvFo = 



(4.12) 



on . However, artificially truncating the initial data 
to compact support causes this integral to be non-zero, 
although small. This leads to the following artificial be- 
havior of Fq at early retarded times. 

This early time behavior of Fq is conveniently de- 
scribed in coordinates {v,l = 1/r) which compactify 
aX I = {]. In these coordinates, the evolution Eq. (4. J) 
becomes 



(2 + 3M/) 



■)df + {1- 3PM)di 



Fo = 0. 



(4.13) 



Consider initial data on X^ with compac t sup port in the 
region Va < v < Vh. Integration of Eq. (4.13) along X^ 
from V = Va to V yields 



diFo 



dvFodv = a{v), 




FIG. 1: Signal -Fb(ii) produced on v = lOOOM by an initial 
compact pulse. The solid curve indicates a fit to a cubic 
polynomial. 



where according to the above argument a{v) should van- 
ish for V > vt- However, our numerical procedure for 
truncating the data leads to a small value of a{vb)- Hence, 
for V > Vh, on any ingoing null hypersurface, Fq has the 
asymptotic behavior Fo{v,l) = a{vb)l + 0{l^). When re- 
expressed in (u, v) double null coordinates, this implies 
Fq has the asymptotic behavior Fq{u,v) = a{vb)u^ + 
0{u^). In numerical simulations, this produces artificial 
cubic polynomial dependence of Fq at early times on X"*" . 
Since this cubic behavior is removed by taking four u- 
derivatives it does not af fect the waveform F4 calculated 
from ^0 using Eq. (3.13). 



We simulated this artificial early time behavior by us- 
ing the advanced time algorithm to propagate an ex- 
tremely narrow initial compact support pulse on X^ of 
the form = [{v — 'yi)('L'2 — v)] for vi < v < V2, and 
Fq = otherwise, where V2 ~ vi — AM. Figure I plots 
the resulting signal Fq on the ingoing null hypersurface 
V = lOOOM. (By looking at the signal on finite v we 
avoid problems associated with spatial infinity i'^.) The 
figure shows that the early time behavior fits a cubic 
polynomial in u. This cubic polynomial dominates at 
later times and masks physical effects such as quasinor- 
mal ringdown. However, Fig. |2| shows that after taking 
four M-derivatives the ringdown is clearly observed and 
there is little trace of the secular behavior. Note the 
drastic difference between the scales of the Figs. |^ and ||. 

In computing waveforms from close approximation 
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FIG. 2: The 4th u-derivative of Fig. 



data, because the initial surface J~ is at finite u the 
artificial early time behavior in _Fo is described by a quar- 
tic rather than cubic polynomial in u but the 4th order 
coefficient is small for reasonably early start times. We 
performed the following experiments to measure the sen- 
sitivity of the Stage II signal to the location of /C~ and 
J~ . In these experiments we used a Stage I run to ob- 
tain close approximation initial data with 77 = 1410 for 
on J~ and measured the signal propagated to I"*" . Ex- 
cept when computational expense would have been pro- 
hibitive, the experiments were carried out with both the 
retarded time and advanced time evolution algorithms 
with virtually equivalent results. 

In the first set of experiments, Kr was placed at 

V = -200M and v = -lOOM, with the initial hyper- 
surface J~ located at it = — 720Af in both cases. For 

V = — lOOM the initial data on J~ implied by the time 
reflection of the Stage I evolution was approximately 5 
orders of magnitude smaller at /C^ than at its peak (at 
■D ~ 0); and for v = — 200M, the size at /C~ was 9 orders 
of magnitude smaller than at the peak. The Stage II 
approximation of smoothing this data to zero on /C^ vio- 
lates the integral condition (4.12) and produces the early 



time polynomial behavior shown in the plots of Fo{u) at 
X+ in Fig. ^. For both locations of K.^ the early time 
behavior is well fit by a 4th order polynomial but the 
size of the polynomial is about 10^ times larger for the 
farther out location of JC^ . 

Figure ^ is a magnified view of Fig. ^ showing how the 
early time quartic polynomial behavior contaminates the 
genuine signal at later times. 

Figure || shows that the two signals agree after taking 
four u-derivatives. 

In the second set of experiments, was located at 
M = -540M, ii = -720M, u = -5090M, and u = -00, 



-800 -600 -400 -200 



FIG. 3: Gomparison of the behavior of Fq for two locations 
of /C". 




-100 



-50 



50 



FIG. 4: Blowup of Fig. ^ in the region of the genuine signal. 



with K, located at v ~ — lOOAf (to accentuate the early 
time polynomial behavior). Figure ^ which plots the 
resulting signal Fo{u) on v = lOOOM, shows that the 
polynomial behavior becomes more pronounced as 
approaches 2^ . 

Figure Q shows how the early time polynomial contam- 
inates the genuine signal at later times, with the effect 
increasing as is located farther in the past. Finally, 
Fig. shows that 9|Fo is unaffected by the artificial 
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FIG. 5: The 4th it-derivative of Fig. ^ shows independence 
of the location of IC~ . The evolution is terminated when late 



time noise begins to contaminate the results. (See Sec. IVE) 




-4500 

-800 -600 -400 -200 



FIG. 6: Comparison of early time behavior for J located at 
four different retarded times. 



FIG. 7: Blowup of Fig. |[ 



FIG. 8: The 4th -u-derivative of Fig |^ showing the indepen- 
dence on the location of . 



E. Late time behavior 



polynomial behavior and indicates that a start time of 
u = — 500M is sufficiently early to obtain the correct 
Stage II retarded waveform on X+ . 

Both sets of experiments confirm that the artificial be- 
havior of Fq at early times has no effect on the accurate 
calculation of the waveform at . 



Both the retarded time and advanced time algorithms 
exhibit exponentially increasing errors at late times. This 
numerical problem arises from the definition of our evo- 
lution variables, which are not designed to deal with late 
time features near the future horizon 7^+. For r w 2M, 
our evolution variables are related to the Israel tetrad 
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variable _Fo by 



4Af2 



Fo, 



(4.14) 



where the aiBne parameter A is exponentially related to 
advanced Bondi time by 

Because Fq is defined in terms of a globally well-behaved 
tetrad, this introduces an exponential dependence in the 
evolution variables which influences the numerical per- 
formance in the following way. 



Near H+, as r ^ 2A/, Eqs. (4.2) and (4.7) reduce to 



-L) 

2M 



dufo = 0, 



(4.15) 



where fo = Fq = IQM'^Fq. This equation admits out- 
going waves of the form fo = /(u) exp(z;/(2M)) in the 
region inside the Schwarzschild potential. These expo- 
nentially growing waves are forbidden in the analytic 
version of our problem by the boundary condition that 
Fo = Fo = on IC~ but they are excited by roundoff 
error. As the computation proceeds, the interval Av be- 
tween IC~ and the peak of the Schwarzschild potential in- 
creases and eventually round-off error in this interval un- 
dergoes sufficient exponential magnification to spoil the 
exterior evolution. There remains an inner region on 7i+ 
where the solution is accurate. For both algorithms, dou- 
bling the precision at which the computation is carried 
out increases the meaningful computation time. For ex- 
ample, the double precision runs shown in Figs. || and 
^ end prematurely due to the onset of numerical noise 
but in quadruple precision both runs can be extended 
through five or six cycles of quasinormal ringdown. 



V. CLOSE APPROXIMATION RETARDED 
WAVEFORMS 

We begin with a review of the close approximation re- 
tarded waveforms from a white hole fission obtained in 
Stage I The close approximation white hole horizon 
perturbation (^41-^-) given by the conformal model is 
a pulse with the following basic properties in terms of 
its dependence on the scale invariant parameter rj (see 
Sec. II A). In the large rj limit {rj > 250) the amplitude 
of the pulse scales quadratically with r] and the width is 
a monotonically decreasing function of rj (the width be- 
comes zero as 77 ^ 00). In the small rj limit (rj < 25) 
the amplitude scales linearly with rj but the shape of the 
pulse is independent of 77. Fig. ^ shows the horizon data 
obtained from various large and mid-size values of 77 (the 
amplitudes have been renormalized so that all pulses have 
the same height). 

The associated retarded white hole fission waveforms 
iFi\x+) have a similar dependence on rj. In the large rj 



T|=7060 : 
Tl=1410 J 1 
Tl=364 I 

Ti=84.3 I I 

I I 



FIG. 9: White hole fission close approximation data F4{u) on 
for 77 = 7060, 1410, 364 and 84.3. 




FIG. 10: Close approximation waveforms F4{u) on T"*" for 
ri = 7060, 1400, 368 and 84.3. 



limit the waveform consists of an extremely sharp ini- 
tial pulse followed by ringdown and power law tail. The 
amplitude of the initial pulse scales quadratically with rj 
and the width is a monotonically decreasing function of 
rj. As 77 is reduced, the initial pulse (signal prior to ring- 
down) becomes wider and obtains more structure. In the 
small 77 limit the amplitude scales linearly with 77 and the 
shape is independent of 77. Fig. |l^ shows a sequence of 
waveforms from large to mid-size values of 77. 
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FIG. 11: Final black hole horizon perturbation after super- 
posing stage I and II, with amplitudes renormalized so that 
the peaks match. The narrow peak corresponds to 77 = 7060 
and r) = 1410. 



We now present the results for the retarded solutions 
from a black hole fusion in the close approximation ob- 
tained by superposition of the advanced black hole results 
(time reversal of the above white hole solutions) and the 
Stage II results. This superposition yields the retarded 
waveform and a modified black hole horizon perturba- 
tion. In the following paragraph, we first describe the 
black hole horizon perturbations. 



Figure 11 shows the perturbation on 7i+ for the re- 
tarded black hole solution. Unlike the Stage I perturba- 
tion of the white hole horizon (see Fig. ^), here the width 
of the perturbation is never smaller than a quasinormal 
period. 

The superposition procedure for obtaining the Stage 
II perturbation of the black hole horizon turns the very 
sharp Stage I pulses associated with large 77 into rela- 
tively broad pulses. The final black hole perturbation 
approaches a unique asymptotic shape in both the large 
and small 77 limits. In Fig. |ll|, rj = 7060 and 77 = 1410 give 
the large 77 shape, while 77 = .00524 and 77 = .00001 give 
the same small 77 shape. Note that the two extremes are 
themselves very similar. Consequently all the retarded 
black hole waveforms obtained from this procedure are 
very similar in shape. 

One of the most striking features of this new hori- 
zon perturbation is the initial ringup, which is later fol- 
lowed by the expected ringdown due to scattering off the 
Schwarzschild potential. The ringup on 7i+ results from 
the ringup of the data on J~ , which in turn results from 
the ringdown of the Stage I waveform. This ringup is 



FIG. 12: Black hole retarded waveforms for close approxima- 
tion data. 



an unavoidable feature of our two stage methodology for 
removing ingoing radiation from the black hole merger. 

Figure ^ shows the waveforms for close approxima- 
tion data obtained from the conformal model. The am- 
plitude of the waveforms scales linearly with 77 and has 
been renormalized so that all waveforms show the same 
quasinormal decay at late times. Note that, aside from 
amplitude, the waveforms have only weak dependence on 
77. This is a great departure from the stage I results that 
the widths of the retarded white hole waveforms depend 
strongly on 77, which produces extremely narrow pulses 
for large 77. 

It should not be surprising that retarded waveforms 
from a black hole merger differ qualitatively from re- 
tarded waveforms from a white hole fission. The fission 
process is directly observable at X+ , whereas the merger 
waveform results indirectly from the black holes through 
the preceding collapse of matter or gravitational wave 
energy that formed them. This also explains why the fis- 
sion waveform should be more sensitive to the parameter 
77 which controls the shape and timescale of the horizon 
data. However, the weakness of the dependence of the 
merger waveform on 77 is surprising. 

The retarded black hole waveforms at X+ also show 
the ringup/ringdown pattern. The ringup is due to 
backscatter of the ringup on the initial slice . Interest- 
ingly, the backscattered signal of an £ = 2 quasinormal 
ringup on J^^ is a ringup on X"*" with the same £ = 2 
time dependence. Figure |l^ shows the ringup phase of 
the waveform obtained by placing JC^ at -D — — lOOM 
and at w = -200Af. (The data for v = -200Af con- 
tains more ringup cycles.) In the figure -F4I1+ is mul- 
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FIG. 13: Ringup of the waveform: F4|2;+ multiplied by 
exp-.08896w/M 



retarded waveform from a binary black hole merger in 
the close approximation. The process of removing ingo- 
ing radiation from the system leads to two major new 
features in the shape of the close approximation wave- 
forms for a head-on collision: (i) an initial quasinormal 
ringup and (ii) weak sensitivity to the parameter -q of the 
model that in some heuristic sense controls the collision 
velocity. Feature (ii) has the potential importance of en- 
abling the design of an efficient template for extracting a 
gravitational wave signal from noise. 

Similar attempts to remove ingoing radiation from 
Cauchy evolutions of close approximation data have not 
yet been made. Such studies would help clarify whether 
the above features are of intrinsic physical origin or an 
artifact of our methodology. 

The present work is part of an ongoing effort to carry 
out a similar two stage characteristic evolution to provide 
the waveforms from coalescing black holes in the nonlin- 
ear regime. Because this is an unexplored area of binary 
black hole physics, this perturbative study of the head- 
on collision provides a preliminary physical check which 
will be useful in extending the work to the nonlinear and 
nonaxisymmetric case, where inspiraling black holes can 
be treated. 



tiplied by exp— .08896w to remove the exponential part 
of the quasinormal time dependence. Note that in the 
V = — 200M case the ringup begins exactly IQQM sooner 
indicating that it is indeed the ringup of the initial data 
that is responsible for the ringup of the waveform on X+. 

VI. DISCUSSION 

We have combined the conformal horizon model with 
a null evolution code to give a new way to calculate the 



Acknowledgments 

We thank Manuela Campanelli for numerous helpful 
discussions. This work has been supported by NSF 
grants PHY 9800731 and PHY 9988663 to the University 
of Pittsburgh. Y. Z. thanks the Albert-Einstein- Institut 
for hospitality. Computer time for this project has been 
provided by the Pittsburgh Supercomputing Center. 



[1] M. Campanelli, R. Gomez, S.Husa, J. Winicour and Y. 

Zlochower, Phys. Rev. D 63 124013 (2001). 
[2] R, H. Price and J. PuUiri, Phys. Rev. Lett. 72, 3297 

(1994). 

[3] Z. Andrado and R, Price, Phys. Rev. D 56, 6336 (1997). 
[4] M. Campanelli, W. Krivan, and C. O. Lousto, Phys. Rev. 

D 58, 024016 (1998). 
[5] M. Campanelli and C. O. Lousto, Phys. Rev. D 58, 

024015 (1998). 

[6] M. Campanelli, C. O. Lousto, J. Baker, G. Khanna, and 
J. PuUin, Phys. Rev. D 58, 084019 (1998). 

[7] R. Glaser, O. Nicasio, R, Price and J. PuUin, Phys. Rev. 
D 59, 044024 (1999). 

[8] W. Krivan and R, Price, Phys. Rev. Lett. 82, 1358 
(1999). 

[9] G. Khanna, J. Baker, R. Glaser, P. Laguna, O. Nicasio, 
H.-P. NoUert. R, Price and J. PuUin, Phys. Rev. Lett. 
83, 3581 (1999). 



[10] M. Campanelh and C. O. Lousto, Phys. Rev. D 59, 
124022 (1999). 

[11] J. Baker, S. Brandt, M. Campanelh, C. O. Lousto, E. Sei- 
del, and R. Takahaslii, Phys. Rev. D 62, 127701 (2000). 

[12] J. Baker, B. Briigmann, M. Campanelli, and C. O. 
Lousto, Class. Quantum Grav. 17, L149 (2000). 

[13] J. Winicour, Prog. Theor. Phys. Supp. 136, 57 (1999). 

[14] R. Gomez, S. Husa, and J. Winicour, Phys. Rev. D 64 
(2001) 024010. 

[15] R. Gomez, S. Husa, L. Lehner, and J. Winicour, Gravi- 
tational Waves from a Fissioning White Hole, in prepar- 
ration. 

[16] N. T. Bishop, R. Gomez, L. Lehner, M. Maharaj, J. Wini- 
cour, Phys. Rev. D 56 6298 (1997). 
[17] R. Gomez, Phys. Rev. D 64 024007, (2001). 
[18] R. Sachs, J. Math. Phys. 3, 908 (1962). 
[19] S. A. Hayward, Class. Quantum Grav. 10, 779 (19993). 
[20] H. Priedrich, Proc. R. Soc. London A 378, 401 (1981). 



14 



[21] A. H. Rcndall, Proc. R. Soc. London A 427, 221 (1990). 
[22] S. A. Toukolsky, Astrophys. J. 185, 635 (1973). 
[23] W. Israel, Phys. Rev. 143, 1016 (1966). 
[24] E. T. Newman and R. Penrose, J. Math. Phys. 3, 566 
(1962). 

[25] L. Lehner, N. T. Bishop, R. Gomez, B. Szilagyi, and J. 

Winicour, Phys. Rev. D 60, 044005 (1999). 
[26] S. Husa and J. Winicour, Phys. Rev. D 60, 084019 

(1999). 

[27] E. T. Newman and T. W. J. Unti, J. Math. Phys. 5, 891 



(1962). 

[28] H. Bondi, M. J. G. van dcr Burg and A. W. K. Metzucr, 

Proc. R. Soc. London A 269 21, (1962). 
[29] R.K. Sachs, Proc. R. Soc. London A 270 103, 1962. 
[30] L. Tamburino and J. Winicour, Phys. Rev. 150, 1039 

(1966). 

[31] B. G. Schmidt, talk at the Tubingen workshop (April, 
2001). 



